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ABSTRACT 


An  analytic  investigation  was  made  to  determine  the 
dynamic  properties  of  a  two-bladed  rigid  fully  gimballed 
helicopter  rotor  incorporating  circulation  control  air¬ 
foils  and  tip  jet  propulsion.  A  time  domain  analysis  was 
developed  which  provided  the  capability  of  using  non¬ 
linear  airfoil  aerodynamics  and  arbitrary  rotor  physical 
characteristics.  The  effects  of  feather  principal  axis 
of  inertia  location,  horizontal  gust  disturbances,  and 
feedback  control  on  rotor  stability  were  assessed. 

Results  of  the  investigation  indicate  that  the  subject 
helicopter  rotors  are  unstable  in  forward  flight  without 
feedback  control.  With  feedback,  the  rotors  are  stable 
and  controllable. 


ADMINISTRATIVE  INFORMATION 

The  research  reported  herein  is  part  of  the  Independent  Research  Program  at  the 
David  W.  Taylor  Naval  Ship  Research  and  Development  Center  (DTNSRDC)  and  during 
fiscal  years  1981  and  1982.  This  research  was  conducted  under  Naval  Material  Command 
Program  Element  61152N,  Task  Area  ZR0230201,  and  Work  Unit  1605-402. 

INTRODUCTION 

The  Aviation  and  Surface  Effects  Department  (ASED)  at  the  David  W.  Taylor  Naval 
Ship  Research  and  Development  Center  (DTNSRDC),  through  its  efforts  to  develop 
advanced  rotorcraft,  has  considered  many  innovative  concepts  which  showed  a  potential 
to  advance  the  rotary  wing  state-of-the-art.  One  such  concept  which  is  believed  to 
have  significant  potential  is  the  fully  gimballed  helicopter  rotor. 

A  gimballed  rotor  is  a  rotor  which  is  fiae  to  pitch  and  roll  at  the  end  of  the 

rotor  shaft  while  it  rotates.  It  can  be  likened  to  a  conventional  rotor  with  a  ball 

joint  replacing  the  rigid  attachment  of  the  rotor  hub  to  the  rotor  shaft  (Figure  1). 

For  a  rotor  to  be  truly  fully  gimballed,  only  forces  can  be  transmitted  between  the 
rotor  and  the  fuselage.  As  a  result,  no  limits  are  placed  on  fuselage  center  of 
gravity  shifts  other  than  those  set  by  fuselage  attitude  considerations.  There  are 
also  no  rotor  control  power  requirements  set  by  center  of  gravity  shifts  because,  in 
steady  flight,  the  resultant  force  on  the  fuselage  always  acts  through  the  rotor 
gimbal.  This  means  that  there  is  a  much  greater  flexibility  in  load  carrying  capa¬ 
bility  for  gimballed  rotors  than  for  other  rotor  types.  It  is  the  major  a^v.ntage 
of  the  concept  and  it  is  particularly  important  for  heavy  lift  helicopters. 


There  are  several  challenges  which  must  be  met  in  order  to  design  a  fully 
gimballed  rotor.  The  first  challenge  is  the  mechanical  design  problem  in  which 
powering  the  rotor  and  providing  collective  and  cyclic  lift  control  must  be 
accomplished  without  passing  torques  between  the  rotor  and  the  fuselage.  Perhaps 
the  most  straightforward  means  of  powering  and  controlling  a  fully  gimballed  rotor 
is  with  a  completely  pneumatic  system  incorporating  tip  drive  propulsion  and  circu¬ 
lation  control.  Tip  drive  systems  have  been  analyzed  extensively,  but  very  few  tip 

1* 

driven  helicopters  have  been  built.  Nichols  gives  a  good  theoretical  analysis  of 

the  efficiency  of  the  pressure  jet  tip  drive  system  in  comparison  with  the  turbine 

2 

powered  shaft  drive,  while  Bossier  and  Harris  compare  various  rotor  drive  systems 
primarily  from  a  weight  standpoint.  Because  there  are  very  few  actual  data  points 
available  for  tip  driven  helicopters,  it  is  difficult  to  draw  hard  and  fast  con¬ 
clusions  as  to  the  superiority  of  drive  system  types.  For  example,  both  Bossier  and 

2  3 

Harris,  and  Head  give  comparative  design  examples,  but  arrive  at  different  con¬ 
clusions  as  to  which  drive  system  is  best.  It  is  clear,  however,  that  tip  drive 
becomes  attractive  for  very  large  helicopters  because  it  is  simple  and  inexpensive 
to  construct. 

Circulation  control  (CC)  is  a  primary  area  of  expertise  in  ASED.  There  are 
4  3  b 

three  major  programs  ’  ’  using  this  technology.  Circulation  control  has  the 
fundamental  advantage  of  providing  a  simple  means  of  modulating  the  magnitude  and 
distribution  of  rotor  lift  from  the  nonrotating  frame  of  a  helicopter.  This  lift 
control  is  accomplished  by  tangentially  blowing  air  through  a  slot  near  the  trailing 
edge  of  the  airfoil.  It  is  a  very  effective  control  and  can  provide  very  high  lift 
coef  f icients. 

The  concept  of  a  two-bladed  fully  gimballed  rotor  incorporating  CC  airfoils 

evolved  from  a  perceived  need  for  a  large  simple  helicopter  with  reduced  ship  board 

stowage  requirements.  Tip  drive  coupled  with  CC  yields  a  mechanically  simple  and 

low  cost  rotor  system.  The  two-bladed  rotor  configuration  enhances  the  design 

simplicity  by  eliminating  the  rotor  folding  requirement  to  meet  stowage  require- 

3 

ments.  In  a  subsequent  design  study.  Head  demonstrated  the  additional  benefit  of 
reduced  gross  weight  for  the  two-bladed  configuration  as  compared  to  higher  numbers 
of  blades.  For  these  reasons,  the  concept  of  a  two-bladed  fully  gimballed  rotor  is 
addressed  in  this  report. 

*A  complete  listing  of  references  is  given  on  page  23. 
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The  second  challenge  which  must  be  met  in  order  to  design  a  fully  gimballed 

rotor  is  to  ensure  that  the  rotor  concept  is  statically  and  dynamically  well  behaved. 

The  best  way  to  meet  this  challenge  is  to  strive  to  gain  a  fundamental  understanding 

of  the  properties  of  this  type  of  rotor.  Two-bladed  gimballed  rotors  have  very  high 

inertias  about  the  rotational  and  flapwise  axes,  but  very  low  inertia  about  the 

feathering  axis.  Thus  the  rotor  has  a  potential  for  instability  especially  if  the 

7 

characteristically  large  pitching  moments  of  CC  airfoils  are  considered.  Chaplin 
addressed  this  configuration  using  a  simplified  analysis  which  assumed  linearized 
aerodynamics,  constant  rotor  speed,  and  zero  or  low  forward  speed.  He  concluded 
that  a  rigid  two-bladed  fully  gimballed  rotor  "can  be  rendered  dynamically  stable  by 
appropriate  mass  balancing." 

This  report  describes  the  methodology  and  results  of  an  effort  to  analyze  the 
dynamic  properties  of  a  rigid  two-bladed  gimballed  rotor  with  circulation  control 
using  a  comprehensive  time  domain  simulation.  No  consideration  is  given  to  fuselage 
coupling  effects,  and  results  should  be  considered  as  equivalent  to  those  obtainable 
in  a  wind  tunnel. 


METHODOLOGY 

Considerable  expertise  has  been  gained  in  time  domain  simulation  as  part  of  the 
CC  rotor  program.  By  using  time  rather  than  frequency  domain  analysis,  it  is 
possible  to  incorporate  more  complete  aerodynamic  and  more  flexible  geometric  repre¬ 
sentations  in  the  simulation. 

A  comprehensive  computer  program  (HE1.STAB) ,  which  simulates  the  blade  and 
vehicle  dynamics  of  articulated  and  hingeless  (A/H)  rotor  helicopters,  was  previously 
developed  by  A3ED  in  support  of  the  CC  rotor  project.  The  general  features  of 
HELSTAB*  are 

1.  Time  history  analysis — to  provide  fully  coupled  degrees  of  freedom,  to 
avoid  the  restrictions  inherent  in  linearized  eigenvalue  solutions,  and  to  allow 
for  a  more  general  vehicle  and  aerodynamic  model. 

2.  Modular  programming — to  maintain  separate  and  distinct  program  modules  for 
ease  of  modification  and  transportability. 


*HELSTAB  has  not  been  formally  documented.  The  primary  contributors  to  its 
development  were  J.B.  Wilkerson,  P.S.  Montana,  and  D.W.  Poe  of  DTNSRDC,  with 
contract  support  from  the  University  of  Maryland  on  specific  tasks,  and  consultative 
support  from  Kaman  Aerospace  in  areas  related  to  the  XH-2/CCR  demonstrator  aircraft. 
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3.  Balanced  approach — to  provide  the  necessary  balance  between  rotor  and 
vehicle  aerodynamic,  dynamic,  and  geometric  descriptions. 

The  gimballed  rotor  is  fundamentally  different  from  A/H  rotors.  This  difference 
stems  from  an  assumption,  common]-  made  in  the  study  of  A/H  rotors,  that  the  rotor 
shaft  angular  velocity  with  respect  to  the  fuselage  is  constant.  Since  the  rotor 
hub  is  rigidly  attached  to  the  rotor  shaft,  the  rotor  blades  may  be  treated  Inde¬ 
pendently  of  each  other  In  solving  for  their  motions.  For  a  gimballed  rotor,  the 
rotor  shaft  velocity  is  not  constant  in  either  magnitude  or  direction.  In  fact,  the 
dynamics  of  gimballed  rotors  are  best  characterized  by  shaft  tilt  angles  rather 
than  by  the  flap,  lag,  and  feather  of  individual  blades.  While  it  may  initially 
appear  that  the  study  of  a  rigid  gimballed  rotor  (i.e.,  no  blade  motion  with  respect 
to  the  rotor  hub)  is  less  general  than  the  study  of  A/H  rotors  with  first  mode  blade 
motion  (i.e.,  rigid  blade  flap,  lag,  and  feather),  the  studies  are  equivalent.  Both 
result  in  first  order  tip  path  plane  motion,  and  both  model  the  n-blade  per  revo¬ 
lution  vibrations  which  affect  the  fupelage. 

The  architecture  of  the  HELSTAB  program  allowed  modifications  for  the  analysis 
of  gimballed  rotors  to  be  incorporated.  The  main  modifications  to  HELSTAB  consisted 
of  t be  replacement  of  program  modules  containing  equations  for  computing:  rotor 
angular  accelerations,  rotor  forces,  and  blade  element  velocities.  Changes  were 
made  to  input  and  output  modules,  and  several  new  modules  were  added  to  the  program 
to  account  for  blade  element  angle-of-attack  variation  due  to  blade  collective  angle 
change  and  for  feedback  of  rotor  angles,  angular  rates,  and  accelerations  to  the 
rotor  control  system.  The  modified  computer  program  (GIMSTAB)  was  exercised  for  a 
variety  of  cases  in  hover  and  forward  flight. 

ROTOR  AERODYNAMIC  REPRESENTATION 

The  HELSTAB  simulation  program  was  designed  to  allow  the  selection  of  the 
levels  of  sophistication  of  rotor  and  vehicle  analyses  via  the  insertion  of 
different  program  modules.  Considerable  effort  has  been  devoted  in  ASED  to  the 
development  of  circulation  control  rotor  (CCR)  aerodynamic  performance  capability. 

The  best  methodology  available  was  incorporated  in  GIMSTAB  when  it  was  developed. 

4 

Wilkerson  et  al  discuss  the  correlation  of  the  rotor  aerodynamic  performance 
methodology  with  full-scale  wind  tunnel  data  obtained  as  part  of  the  CCR  flight 
demonstrator  test  program.  In  brief,  the  aerodynamic  model  is  based  on  strip  theory 


4 


with  Glauert  inflow.  All  CC  airfoil  characteristics  (e.g.,  chord,  thickness,  slot 
height)  could  be  varied  radially.  The  basic  airfoil  characteristics  were  input  as  a 
series  of  data  tables  for  the  specific  airfoils  used.  Corrections  were  made  for  tip 
losses,  and  compressibility  and  Reynolds'  effects.  The  model  was  used  for  calcu¬ 
lating  the  radial  distribution  of  blowing  jet,  momentum  coefficient,  and  weight  flow 
in  addition  to  the  lift,  drag,  and  pitching  moments.  The  blade  lift,  drag,  and 
pitching  moments  were  provided  as  inputs  to  the  rotor  dynamics  portion  of  the 
GIMSTAB  program. 


EQUATIONS  OF  MOTION 

The  equations  of  motion  for  rigid  gimballed  rotors  were  derived  in  a  relatively 
conventional  fashion  based  on  Newton's  Second  Law  and  closely  following  the  method 

g 

and  notation  of  Etkin.  The  analytic  rotor  model  selected  was  very  general  and 
provided  a  maximum  flexibility  in  the  description  of  the  rotor.  Derivations  were 
made  with  the  following  limitations  placed  on  the  model  and  program  architecture. 

1.  The  rotor  was  rigid.  The  motion  of  the  rotor's  center  of  gravity  and  of 
each  blade  is  described  explicitly  by  the  motion  of  the  gimbal  in  inertial  space 
plus  the  kinematics  of  the  rotor  about  the  gimbal. 

2.  The  motion  of  the  gimbal  was  prescribed  by  a  global  vehicle  dynamics 
solution  external  to  the  rotor  dynamics  model. 

3.  Total  aerodynamic  forces  and  moments  were  computed  and  supplied  to  the  rotor 
dynamics  model  by  a  separate  aerodynamics  model. 

4.  Tip  jet  thrust  was  also  supplied  to  the  dynamics  model. 

5.  Gimbal  reaction  forces  were  computed  by  the  dynamics  model. 

6.  Gimbal  restraint  moments  were  accommodated  within  the  dynamics  model  to 
provide  for  frictional  losses  and  the  potential  need  for  damping  and/or  centering 
springs  as  rotor  motion  limiters.  These  provisions  were  in  the  nature  of  damping 
and  friction  coefficients  and  spring  rates  coupled  to  angular  rates  and  deflections 
to  produce  torques  about  the  gimbal.  The  nature  of  the  structural  limitations 
imposed  on  the  model  were  in  keeping  with  the  modular  nature  of  the  GIMSTAB 
program. 

The  equations  of  motion  of  the  rotor  were  derived  in  a  reference  frame  (Figure 
2)  attached  to  the  rotor  located  at  its  center  of  gravity  (c.g.).  In  this  reference 
frame,  rotor  inertias  are  invariant  and  rotor  weight  does  not  contribute  a  moment. 
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A  draw  back  of  the  c.g.  centered  frame  is  the  appearance  of  gimbal  reaction  forces 
in  the  moment  equations.  However,  the  kinematic  relationship  of  linear  motion  to 
angular  motion  (i.e.,  rigid  rotor  assumption)  allows  the  force  equations  to  be  used 
to  determine  gimbal  reaction  rather  than  rotor  c.g.  linear  motion.  Thus,  the  six 
fundamental  equations  of  motion  were  reduced  to  three  moment  equations.  As  a 
consequence  of  the  time  history  approach,  these  three  equations  were  solved 
explicitly  for  three  Euler  angular  accelerations.  (Euler  angles  used  in  the 
equation  of  motion  are  shown  in  Figure  3.)  A  complete  description  of  the  rotor 
motion  was  then  available  by  numerically  integrating*  the  angular  accelerations. 

The  details  of  the  derivation  of  the  equations  of  motion  are  as  follows.  The 
symbols  are  defined  in  the  separate  notation  section;  however,  the  meaning  of  most 
symbols  can  be  determined  from  the  equations  and  from  Figure  4. 

Newton’s  Second  Law  yields  six  equations 


Rate  of  Change 
of  Angular 
.Momentum 

♦ 

h 

loo  +  (0  x  1(0 


Rate  of  Change 
of  Linear 
.  Momentum 


Sum  of  the 
Applied  External 
.Moments 

G 

G 


Sum  of  the 
Applied  External 
.Forces 


(1) 


F 


which  reduces  to 


mr 


F 


(2) 


for  constant  rotor  mass.  The  total  external  moment  applied  to  the  rotor  consists  of 
*The  integrations  were  carried  out* using  a  four-step  Runge-Kutta  (R-K)  formu- 

Q 

lation.  This  technique  was  compared  to,  a  six-step  R-K  and  variable  time  step 
methods  and  was  found  to  have  adequate  numerical  accuracy  and  the  highest  time 
efficiency. 
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Iw  +  co  x  Ico  **  Ma  -b  x  (ir’-g)  m  +  (6) 

IL  ir  this  final  set  of  Equations  (6)  which  must  be  solved  for  the  Euler  angle 
••  ••  *•  • 

accelerations  (<j>,0,^).  These  angular  accelerations  appear  in  only  two  terms,  oj  and 
if.  The  solut' ion  procedure  is  shown  symbolically  in  Equations  (7 )  —  (10)  in  which  the 
previous  vector  notation  is  replaced  with  matrix  notation  and  where 
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IGJ  -  [A] 


r  -  [a] 


+  {D} 


+  {d) 
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$ 

+  {d}  +  [Oil  LI]  {w}  -  {Ma}  -  [b]  [a]  J  0 

ij 


+  (d}  -  {g }  m  +  (MR}  +  {MTJ}  (7) 


[A]  +  [b]  [a]  m 
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JMS}  -  [b]  Ld)  -  {g}  m  +  (MR>  -  (P)  -  [co]  [I]  (co)  4-  (8) 


[K] 


(C) 


(9) 


[K]"1  {C}  (10) 

It  la  worth  noting  here  that  the  solution  was  a  complex  undertaking.  The  main 
difficulty  was  associated  with  the  intermediate  frames  of  reference  between  the 
glmbal  and  the  rotor  in  which  the  Euler  angles  (<t>,0,i|>)  were  defined.  This  inherent 
cross  coupling  Introduced  a  myriad  of  terms  consisting  of  products  of  direction 
cosines.  No  small  angle  assumptions  were  made  to  simplify  the  equations.  The 
solution  was  accomplished  with  the  use  of  the  Project  MAC's  Symbolic  Manipulation 


l 


8 


System"^  (MACSYMA)  available  on  the  Navy  Laboratory  Computer  Network.  The  main 
benefit  of  MACSYMA  in  this  application  was  its  ability  to  symbolically  manipulate 
matrices.  It  was  also  possible  to  store  results  and  to  automatically  convert 
equations  to  FORTRAN  coding  on  the.  system.  With  intermediate  steps,  the  FORTRAN 
statements  were  transmitted  to  DTNSRDC's  computers  for  incorporation  in  the  GIMSTAB 
program. 


DISCUSSION 

The  objective  of  this  program  was  to  analyze  the  dynamic  properties  of  two- 
bladed  rigid  glmballed  rotors  incorporating  CC  airfoils.  Time  domain  analysis  was 
selected  to  accomplish  the  objective  to  allow  the  use  of  nonlinear  aerodynamics  and 
to  provide  a  flexible  and  complete  geometric  model  of  the  rotor.  There  were  three 
steps  taken  to  accomplish  the  objective:  (1)  selection  of  a  rotor  physical  model, 
(2)  determination  of  hover  characteristics,  and  (3)  simulation  of  forward  flight 
characteristics. 

ROTOR  MODEL 

At  the  time  the  the  glmballed  rotor  project  was  executed,  the  best  defined  CC 
rotor  was  the  one  employed  by  the  XH-2/CCR  flight  demonstrator  aircraft.  The 
airfoil  characteristics  were  well  known  and  data  tables  were  in  existence  for  use 
in  rotor  performance  simulations.  In  addition,  the  actual  physical  characteristics 
of  the  flight  demonstrator’s  rotor  had  been  measured  and  were  available.*  The 
physical  characteristics  of  the  XH-2/CCR  rotor  blades  and  hub  were  modified  to 
reflect  the  change  from  four  blades  to  two  blades;  an  emphasis  was  placed  on 
preserving,  as  much  as  possible,  the  combined  hub  and  blade  characteristics.  For 
instance,  total  hub  and  blade  weight  and  rotational  inertia,  and  solidity  were  held 
constant.  Individual  blade  characteristics  were  scaled  as  follows: 

1.  Radial  location  of  the  blade  center  of  gravity  was  held  constant. 

2.  Blade  thickness  ratio  distribution  was  held  constant. 

3.  Blade  weight  was  doubled. 

4.  Blade  chord  was  doubled. 

5.  Blade  inertias  were  scaled  consistent  with  items  1  through  4. 


*These  characteristics  were  provided  by  Kaman  Aerospace  Corporation. 
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6.  Blade  airfoil  selection  was  held  constant. 

7.  All  blade  nondimens ional  CC  parameters  were  held  constant. 

A  description  of  the  rotor  model  is  given  in  Tables  1  and  2. 

HOVER  ANALYSIS 

The  hover  analysis  established  the  baseline  characteristics  of  the  gimballed 
rotor  model.  The  central  feature  of  these  characteristics  was  the  development  of  a 
stability*  boundary  in  regard  to  the  location  of  the  rotor's  principal  axis  of 
inertia.  As  noted  in  the  introduction,  a  two-bladed  gimballed  rotor  has  very  high 
inertias  about  two  axes,  flap  arid  azimuth,  and  very  low  inertia  about  one  axis, 
feather.  Chaplin  identified  the  feather  principal  axis  of  inertia  location 
relative  to  the  blade  aerodynamic  center  as  a  key  parameter  affecting  rotor  dynamic 
stability  for  two-bladed  gimballed  rotors.  Feather  principal  axis  location  is 
important  whether  CC  or  conventional  airfoils  are  employed;  however,  the  use  of  CC 
airfoils  complicates  the  relationship  of  principal  axis  of  inertia  location  to 
stability. 

A  salient  feature  of  circulation  control  airfoils  is  the  existence  of  two 
aerodynamic  centers.  Conventional  airfoils  employed  on  helicopters  have  one  aero¬ 
dynamic  center  (a.c.)  located  near  the  quarter  chord.  Since  angle  of  attack 
determines  lift  (for  a  given  velocity  and  air  density),  the  location  of  the  principal 
axis  (center  of  gravity  for  sections)  with  respect  to  the  a.c.  determines  whether 
the  response  to  an  angle-of-attack  disturbance  is  stable  or  unstable.  For  stability, 
the  principal  axis  should  be  ahead  (upstream)  of  the  a.c.;  thus,  a  positive  lift 
disturbance  would  yield  a  nose  down  moment  tending  to  reduce  angle  of  attack  and 
lift. 

Circulation  control  airfoils  have  both  angle  of  attack  and  CC  lift.  The  CC 
lift  is  relatively  independent  of  angle  of  attack,  but  is  dependent  upon  the 
characteristics  of  the  trailing  edge  blowing.  The  a.c.  for  CC  lift  is  located  near 
the  half  chord.  When  the  complex  flow  over  a  rotor  in  forward  flight  is  considered, 
it  is  apparent  that  the  combined  lift  on  the  individual  airfoil  sections  along  a 
blade  may  act  at  the  quarter  chord,  half  chord,  or  anywhere  in  between. 


*The  term  "stability"  is  used  in  this  report  in  several  ways,  some  of  which  do 
not  precisely  match  traditional  interpretations.  Definitions  of  "stability"  used  in 
this  report  are  given  on  pages  11  and  14. 
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In  order  to  classify  feather  principal  axis  locations,  the  0.7  rotor  radius 
station  was  selected  as  a  reference  for  axis  location.  The  location  is  specified  as 
a  fraction  of  chord  length  from  the  leading  edge.  The  0.7  radius  station  is  a 
standard  aerodynamic  reference  and  is  considered  to  be  typical  of  rotor  aerodynamics 
based  on  strip  theory  combined  with  nominal  tip  losses. 

Part  of  the  definition  of  rotor  characteristics  was  the  generation  of  a  rotor 
performance  map  in  hover.  The  map  covers  a  range  of  parameters  representative  of 
the  XH-2/CCR  aircraft.  The  parameters  include:  thrust  from  0  to  13,000  lb 
(57,82.9  N),  blade  duct  pressure  ratios  from  1.0  to  1.8,  and  blade  collective  feather 
angles  from  -6  to  +10  deg  (-0.10  to  0.017  rad).  The  map  was  generated  for  a 
principal  axis  of  inertia  location  at  30-percent  chord  at  the  0.7  radius  station. 

An  assumption  made  during  the  hover  performance  analysis  was  constant  average  tip 
speed  of  615  ft/sec  which  was  accomplished  through  the  use  of  constant  tip  jet 
thrust.  Tip  jet  thrust  depended  upon  blade  collective  angle,  duct  pressure  ratio, 
and  rotor  thrust. 

Stability  in  hover  was  defined  as  the  ability  to  sustain  a  vertical  rotor  axis 
position  during  time  integration  without  simulated  disturbances.  While  this  was  an 
unsophisticated  test,  it  did  succeed  in  establishing  a  boundary  which  was  later 
correlated  by  the  method  of  Chaplin.^  Chaplin's  method  was  based  on  constant 
instantaneous  rotor  tip  speed.  Hence,  his  stability  boundary  is  not  a  function  of 
the  distribution  of  lift  between  angle  of  attack  and  circulation  control.  The 
present  results  confirm  this  assumption  by  showing  only  a  small  sensitivity  to  angle 
of  attack  versus  CC  lift  distribution.  There  is  a  need  to  investigate  the  possibil¬ 
ity  of  stable  precessions!  modes  in  hover  which  is  not  addressed  here. 

FORWARD  FLIGHT 

The  analysis  of  the  motions  of  helicopters  differs  from  that  of  fixed  wing 
aircraft  in  the  ability  to  simplify  the  problem  by  separating  modes.  For  fixed-wing 
aircraft,  it  is  possible  to  consider  the  lateral  and  longitudinal  motions  separately 
because  of  the  existence  of  left  to  right  geometric  and  flow  symmetries.  Helicopters 
possess  the  same  type  of  geometric  symmetry  in  the  fuselage  and  have  geometric 
rotational  symmetry  of  the  rotor  blades  about  the  rotor  hub.  However,  the  flow  over 
a  helicopter  is  not.  symmetric  in  either  hover  or  forward  flight.  Rotational  velocity 
components  exist  in  the  rotor  downwash  which  affect  the  sides  of  the  fuselage 
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differently.  In  forward  flight,  the  velocity  of  each  blade  element  is  composed 
mainly  of  forward  speed  and  rotational  speed;  hence  it  is  a  function  of  the  azimuth 
position  of  the  rotor  blade.  The  blades  of  most  rotors  are  treated  individually. 

In  steady  state  forward  flight,  they  experience  a  once  per  revolution  cyclic  aero¬ 
dynamic  perturbation.  The  blades  of  a  rigid  gimballed  rotor  are  treated  collectively 
because  their  motions  are  kinematically  related.  As  a  result,  the  whole  rotor, 
including  each  blade,  responds  to  a  cyclic  aerodynamic  perturbation  occurring  at  the 
rotor's  rotational  speed  times  the  number  of  blades,  n,  per  revolution.  Hence,  in 
steady  state  forward  flight,  the  dynamics  of  a  gimballed  rotor  can  be  expected  to 
repeat  every  1/n  revolution  (i.e.,  at  the  same  frequency  as  the  forcing  function). 

The  normal  means  of  determining  whether  or  not  an  articulated  rotor  has 
reached  blade  dynamic  equilibrium  in  the  time  domain  is  to  integrate  the  rotor 
motions  for  one  complete  revolution  and  compare  the  starting  and  ending  blade 
angles  of  flap,  lag,  and  feather.  If  the  angles  repeat  within  specified  tolerances, 
the  dynamics  are  defined  as  being  in  equilibrium.  This  method  is  possible  because 
of  the  tacit  assumption  of  constant  rotor  speed.  This  assumption  makes  time  and 
rotor  azimuth  angle  equivalent.  The  individual  rotor  blades  also  experience  a  very 
strong  centrifugal  force  which  tends  to  quickly  stabilize  the  motions. 

For  gimballed  rotors,  the  dynamic  angles  are  rotor  roll,  pitch,  and  azimuth. 
Since  rotor  speed  (i.e.,  azimuth  angle  rate)  is  not  constant,  the  normal  means  for 
determining  dynamic  equilibrium  is  not  possible  with  a  constant  time  step  integration 
technique.  The  practical  alternative  to  a  sophisticated  variable  time  step  inte¬ 
gration  was  an  average  rotor  speed  controller.  This  approach  is  also  realistic  be¬ 
cause  of  the  aerodynamic,  and  structural  design  optimization  of  all  rotors  for  a  given 
rotor  speed. 

A  twofold  approach  was  used  to  achieve  constant  average  rotor  speed.  The 
azimuth  position  at  the  end  of  the  integration  interval  (time  for  one-half 
revolution  of  a  two-bladed  rotor  at  the  desired  average  rotor  speed)  was  compared 
to  the  starting  position.  If  the  rotor  had  not  completed  or  had  exceeded  a  half 
revolution,  the  average  tip  jet  thrust  was  adjusted  to  correct  the  error.  During 
the  time  integration,  the  rotor  azimuthal  acceleration  was  monitored;  and  the 
average  tip  jet  thrust  was  corrected  (with  a  maximum  of  50%  authority)  tc  attempt  to 
achieve  zero  acceleration.  This  unsophisticated  method  proved  to  be  fairly  success¬ 
ful  when  the  rotor  was  stable  or  mildly  unstable  (i.e.,  near  equilibrium). 
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In  spite  of  the  use  of  rotor  speed  control,  it  took  lengthy  simulations  and 
many  iterations  of  cyclic  blowing  to  reach  equilibrium  for  a  specified  rotor 
collective  angle  and  thrust  value.  The  strong  stabilizing  influence  of  centrifugal 
force  is  not  present  for  gimballed  rotors.  By  nature,  gimballed  rotors  are  different 
from  articulated  rotors.  Rotor  dynamic  equilibrium  is  a  much  more  restrictive 
condition  for  the  gimballed  rotor,  because  it  requires  that  the  net  moments  acting 
on  the  rotor  be  zero  in  addition  to  the  other  common  requirements. 

FEEDBACK  CONTROL 

As  a  result  of  experience  using  the  simulation  program,  it  was  obvious  that  a 
means  of  reaching  dynamic  equilibrium  in  forward  flight  was  needed.  After  several 
unsuccessful  trials  using  rotor  tilt  angle  (pitch  and  roll)  feedback,  it  was  found 

that  feeding  back  of  rotor  teeter  rate  was  a  highly  effective  means  of  quickly 

i* 

reaching  dynamic  equilibrium.  The  algorithm  developed  for  this  type  of  feedback 

wa3  very  general  in  nature  and  was  capable  of  generating  feedback  signals  to  blade 

TT 

duct  pressure  for  teeter  and  feather  angles  and  their  first  and  second  derivatives. 

ftt 

The  relationship  used  to  generate  a  control  pressure  signal  was 

Ap  -  kg  Ag*  +  k£  Ag  +  kg  Ag*  +  kg  A0*  +  A0*  +  k^  A0* 


where  each  of  the  A-quantlties  (with  the  exception  of  pressure)  was  the  difference 
between  actual  (simulated)  and  commanded  angles,  rates,  or  accelerations.  The  k- 
quantities  were  simple  gain  constants.  Feedback  was  implemented  by  modifying  the 
normal  blade  duct  pressure  relation. 


duct 


collective 


+  A  sin  +  B  cos  ^ 


The  use  of  teeter  rate  feedback  was  an  extension  of  a  suggestion  made  by 
H.R.  Chaplin,  in  informal  communications,  that  the  use  of  teeter  angle  feedback 
might  be  effective  in  stabilizing  the  rotor. 

rThe  teeter  and  feather  angles  are  Euler  angles  which  are  indicated  in 
Figure  5. 

ftf 

Because  the  control  pressure  signal  is  a  function  of  feather  and  teeter 
angles,  which  are  harmonic  functions  of  azimuth  angle,  the  control  pressure  signal 
is  also  a  harmonic  function  of  azimuth  angle. 
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using  the  new  cyclic  pressure  coefficient  definitions 


A  ■  A^  +  Ap  sin  i|>* 

B  ■  +  Ap  cos  ip* 

where  A^  and  B^  were  the  previous  (commanded)  values.  Substituting  the  new 
definition  results  in 


P 


duct 


^collective  +  A1  s*n  ^  +  B1  C0S  ^  +  AP(sin  sin  'P  +  cos  cos 


or 


’duct  ~  Collective  +  A1  sin  ^  +  Bi  cos  ^  +  AP 


(for  small  values  of  B  and  0) . 


STABILITY 

Several  tests  were  used  to  establish  rotor  stability  characteristics.  The  first 
test  was  whether  or  not  rotor  dynamic  equilibrium  could  be  reached  without  feedback. 
This  was  a  very  cumbersome  process  since  the  required  cyclic  pressure  and  rotor 
collective  angle  combination  for  a  desired  rotor  tilt  angle  was  not  known,  and  the 
rotor  had  to  be  very  close  to  the  equilibrium  dynamics  to  begin  to  stabilize.  The 
second  test  was  to  disturb  the  rotor  from  equilibrium  with  a  pulsed  horizontal  gust 
of  one-half  revolution  duration.  The  pulse  test  was  used  to  determine  if  the  rotor 
possessed  any  inherent  damping.  The  third  test  was  to  disturb  the  rotor  from 
equilibrium  with  a  step  horizontal  gust.  The  intent  of  this  test  was  to  determine 
if  the  rotor's  response  was  predictable  and  if  the  rotor  would  transit  to  a  new 
equilibrium  position.  All  tests  were  performed  both  with  and  without  feedback. 


RESULTS 

HOVER  ANALYSIS 

The  hover  analysis  yielded  two  products,  a  rotor  performance  map  and  a  principal 
axis  stability  boundary.  The  rotor  performance  map  (Figure  6)  was  generated  for  a 
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range  of  parameters  representative  of  the  XH-2/CCR  aircraft:  thrust  from  0  to 
13,000  lb  (57,824  N),  blade  duct  pressure  ratios  from  1.0  to  1.8,  and  blade 
collective  angles  from  -6  to  +10  deg  (-0.105  to  0.175  rad).  The  rotor  feather 
principal  axis  was  located  at  30%  chord  at  the  0.7  radius  station.  The  map  shows 
that  a  desired  rotor  thrust  can  be  achieved  with  a  variety  of  combinations  of  blade 
duct  pressure  ratio  and  blade  collective  angle.  It  also  shows  that  tip  jet  thrust 
can  be  minimized  by  selecting  the  proper  combination  of  the  same  variables  for  a 
given  rotor  thrust  setting. 

The  rotor  principal  axis  stability  boundary,  as  discussed  previously,  is  a 
boundary  describing  conditions  for  which  rotor  dynamic  equilibrium  in  hover  at  a 
thrust  of  11,000  lb  (48,928  K)  was  possible.  The  boundary  (Figure  7)  proved  to  be 
insensitive  to  the  distribution  of  lift  between  angle  of  attack  and  CC  until  a  blade 
duct  pressure  ratio  of  about  1.7  was  exceeded.  This  value  of  duct  pressure  ratio 
corresponds  approximately  with  the  region  for  which  angle-of-attack  lift  changes 
from  positive  to  negative.  It  is  not  unreasonable  to  expect  a  change  in  stability 
characteristics  near  this  condition. 

FORWARD  FLIGHT 

Forward  flight  analysis  proved  to  be  more  complicated  than  originally  en¬ 
visioned.  Achieving  dynamic  equilibrium,  with  the  attendant  requirement  of  constant 
average  rotor  speed,  was  the  most  significant  problem.  Once  the  difficulty  reaching 
dynamic  equilibrium  was  fully  appreciated  the  other  areas  of  Investigation  were 
developed  as  logical  extensions  of  the  effort. 

The  forward  flight  analysis  was  limited  to  a  forward  speed  of  50  ft/sec 
(15.24  m/s)  and  a  thrust  level  of  12,000  lb  (53,376  N) . 

Open  Loop  Dynamic  Equilibrium 

The  conditions  for  which  rotor  dynamic  equilibrium  was  desired  were  -5  deg 
(-0.087  rad)  rotor  forward  tilt  and  0.0  deg  (0  rad)  lateral  tilt  with  the  feather 
principal  axis  located  at  x/c  -  0.2.  This  feather  principal  axis  location  was 
selected  because  of  stability  considerations  which  indicated  that  the  center  of  lift 
on  the  rotor  blade  should  always  be  located  aft  of  the  principal  axis.  Blade  duct 
pressure  was  selected  as  1.8  to  allow  for  the  potential  of  maximum  cyclic  blowing. 
Equilibrium  dynamics  were  achieved  using  an  exhaustive  trial  and  error  method. 


The  required  cyclic  blowing  pressure  ratios  were  (A,B)  =  (-0.0327,  0.0384)  which 
were  very  nominal  values;  the  required  collective  pitch  angle  was  -4.6  deg 
(-0.080  rad).  Figure  8  shows  about  14  consecutive  rotor  revolutions.  A  small  amount 
of  drift  was  still  present  in  the  two  per  revolution  rotor  "wobble"*;  however,  the 
conditions  were  sufficiently  stable  for  equilibrium  to  be  assumed. 

Effect  of  Principal  Axis  Location 

The  effect  of  principal  axis  location  was  evaluated  for  two  cases:  100%  angle- 
of-attack  lift  and  about  100%  CC  lift.  In  both  cases,  the  feather  principal  axis 
location  was  varied  from  20%  to  40%  of  the  chord  at  the  0.7  radius  station.  Non¬ 
equilibrium  initial  conditions  were  selected  to  amplify  the  dynamic  tendencies  of 
the  rotor.  The  initial  conditions  were:  rotor  tilt  angles  equal  to  zero;  rotor 
azimuth  rate,  27.9545  rad/s;  all  other  rates  and  accelerations  equal  to  zero.  No 
cyclic  blowing  was  used.  For  the  100%  angle-of -attack  lift  cases,  the  collective 
angle  was  7.5  deg  (0.131  rad)  and  the  collective  pressure  ratio  was  1.0;  and  for  the 
100%  CC  lift  cases,  they  were  -4.6  deg  (-0.080  rad)  and  1.8,  respectively. 

The  results  are  shown  in  Figures  9  and  10  as  plots  of  rotor  tilt  angles,  $ 
versus  G.  These  plots  show  that  the  rate  of  departure  of  the  rotor  from  the  non¬ 
equilibrium  starting  condition  and  the  frequency  and  amplitude  of  the  rotor  motions 
are  functions  of  the  feather  principal  axis  location  and  distribution  of  lift  be¬ 
tween  CC  and  angle  of  attack,  specifically: 

1.  The  100%  CC  lift  cases  depart  from  the  initial  condition  at  less  than 
one-half  the  rate  of  the  equivalent  100%  angle-of-attack  lift  cases.  This  result 
reflects  the  postulate  that  CC  lift  is  insensitive  to  angle  of  attack.  (Angle  of 
attack  and  local  velocity  vary  with  azimuth  position  in  forward  flight.) 

2.  The  minimum  departure  rates  for  both  the  CC  and  angle-of-attack  lift  cases 
occur  for  a  feather  principal  axis  location  of  30%  chord.  This  result  indicates 
that  there  is  an  optimum  principal  axis  location  for  each  CC  rotcr  design  and  that, 
contrary  to  the  original  belief,  a  too  far  forward  principal  axis  location  can 
reduce  rather  than  enhance  stability. 


*"Wobble"  refers  to  the  conical  motion  that  the  axis  (Figure  2)  attached  to 

the  rotor  hub  would  scribe.  Refer  also  to  Figure  11a  which  shows  rotor  angular 
motions  in  the  more  traditional  form  as  plots  versus  azimuth  angle. 
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3.  The.  minimum  amplitude  of  the  tilt  angle  oscillations  occur  for  a  feather 
principal  axis  location  of  30%  chord. 

4.  The  frequency  of  the  tilt  angle  oscillations  changes  from  two  to  one  per 
revolution  as  the  principal  axis  location  moves  from  25%  to  40%  chord.  The  effect 
is  shown  in  Figure  11. 

FEEDBACK  CONTROL 

All  six  types  cf  feedback  provided  in  the  duct  pressure  feedback  relationship 
were  evaluated  individually.  Initially,  the  purpose  of  feedback  was  to  accelerate 
the  convergence  of  rotor  motion  to  dynamic  equilibrium.  The  test  used  to  evaluate 
the  effect  of  feedback  was  to  start  the  rotor  in  an  unstable  condition  and  monitor 
its  progress  toward  equilibrium.  All  tests  were  performed  using  the  same  starting 
conditions — zero  rotor  pitch  and  roll  angles,  27.9545  rad/s  azimuth  rate,  50  ft/sec 
(15.24  m/s)  flight  velocity,  zero  cyclic  and  1.8  collective  pressure  ratios,  -4.6 
deg  (-0.0803  rad)  collective  angle,  and  20%  chord  feather  principal  axis  location. 

The  feedback  control  signals  were  generated  by  comparing  the  rotor's  computed 
angles,  rates,  or  accelerations  with  those  equivalent  to  the  average  open  loop 
equilibrium  condition  of  -5  deg  (-0.0873  rad)  pitch  and  0  deg  (0  rad)  roll.  The 
results  of  the  feedback  trials  (Figures  12  and  13)  were  compared  with  the  open  loop 
results  for  the  same  starting  conditions  (Figure  9a)  to  obtain  a  qualitative  assess¬ 
ment  of  the  value  of  feedback. 

Teeter  (Flap)  Feedback 

The  effects  of  the  various  types  of  teeter  feedback  are  shown  in  Figure  12. 

Both  positive  and  negative  gain  constants,  kg,  k»,  or  kg,  were  used.  The  results 
are  summarized  as  follows: 

1.  Negative  teeter  angle  feedback  is  slightly  stabilizing  in  that  the  resulting 
complex  dynamics  appear  to  be  centered  near  the  desired  roior  tilt  angles.  The 
dynamics  contain  a  two  per  revolution  oscillation  with  variable  amplitude  super¬ 
imposed  on  a  longer  period  oscillation  whose  frequency  is  a  function  of  gain 

constant,  k„ . 

P 

2.  Positive  teeter  angle  feedback  is  destabilizing. 

3.  Positive  teeter  rate  feedback  is  very  stabilizing.  For  k^  ■  400,  the 
rotor  reached  a  near  equilibrium  dynamic,  condition  within  about  1.5  seconds.  The 
rotor  oscillations  were  two  per  revolution  with  constant  amplitude. 
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4.  Negative  teeter  rate  feedback  is  destabilizing. 

5.  Negative  teeter  acceleration  feedback  is  slightly  stabilizing  in  the  sense 
that  the  departure  rate  from  the  starting  condition  is  reduced  from  the  open  loop 
case.  The  character  of  the  oscillations  is  two  per  revolution  with  near  constant 
amplitude.  Large  values  of  teeter  acceleration  feedback  gain,  k^,  change  the 
character  of  the  oscillations  to  one  per  revolution  with  increasing  amplitude  and 
hence,  have  a  destabilizing  influence. 

6.  Positive  teeter  acceleration  feedback  is  destabilizing. 

Feather  Feedback 

The  effects  of  various  types  of  feather  feedback  are  shown  in  Figure  13.  The 
results  are  summarized  as  follows: 

1.  Negative  feather  angle  feedback  is  slightly  stabilizing  in  that  it  reduces 
the  rate  of  departure  from  starting  conditions  as  compared  to  the  open  loop  case. 

2.  Positive  feather  angle  feedback  is  destabilizing. 

3.  Negative  feather  rate  feedback  is  stabilizing;  there  appears  to  be  an 
optimum  gain  value  near,  k^  =  -100.  Larger  gain  values  seem  to  excite  lower  fre¬ 
quency  oscillation.  Characteristics  for  k^  *  -100  include  the  expected  two  per 
revolution  oscillations  and  relatively  steady  dynamics. 

4.  Positive  feather  rate  feedback  is  destabilizing. 

5.  Negative  feather  acceleration  feedback  did  not  appear  to  be  stabilizing  for 
the  length  of  the  simulation. 

6.  Positive  feather  acceleration  feedback  is  destabilizing. 

Effect  of  Lift  Type 

Several  trials  were  made  to  assess  the  effect  of  the  type  of  lift,  angle  of 
attack,  or  CC  on  rotor  dynamic  characteristics  with  teeter  rate  feedback,  k£  =  400. 
The  results  of  these  trials,  which  had  the  same  initial  conditions  as  above,  are 
shown  in  Figure  14.  In  general,  the  distribution  of  lift  between  types  did  not 
affect  the  ability  of  the  rotor  to  reach  dynamic  equilibrium  or  the  speed  at  which 
equilibrium  was  approached.  The  type  of  lift  distribution  did  influence  the 
location  of  equilibrium  dynamics  slightly.  This  is  shown  graphically  in  Figure  15 
as  plots  of  mean  rotor  equilibrium  pitch  and  roll  tilt  angles  versus  blade  collective 
pressure  ratio.  There  is  a  previously  unmentioned  feature  of  the  simulations 
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pneumatic  algorithm  which  should  be  noted.  The  pneumatic  algorithm  was  designed 
to  prevent  the  evaluation  of  rotor  aerodynamics  for  duct  pressure  ratios  less  than 
1.0,  a  condition  which  corresponds  to  suction  rather  than  blowing  at  the  CC  slot. 
While  suction  can  be  a  very  powerful  boundary  layer  control,  it  is  not  the  phenomena 
which  CC  airfoils  were  designed  to  exploit.  No  data  are  available  for  CC  airfoils 
with  suction  instead  of  blowing.  And,  in  reality,  it  would  be  difficult  to  design 
a  pneumatic  system  to  accommodate  mass  flow  in  either  direction.  For  these  reasons, 
whenever  a  duct  pressure  ratio  less  than  1.0  was  calculated,  the  duct  pressure  ratio 
was  assigned  a  value  of  1.0. 

For  the  trials  under  consideration,  when  the  collective  pressure  ratio  was  1.0 
(i.e.,  the  rotor  only  had  angle-of-attack  lift),  the  feedback  control  algorithm  was 
only  allowed  to  generate  a  pressure  signal  when  the  change  in  pressure  was  positive. 
This  is  the  probable  cause  for  the  shift  in  mean  rotor  roll  tilt  angle  as  the 
pressure  ratio  approached  1.0  which  was  evident  in  Figure  15.  The  ability  of  the 
rotor  to  stabilize  with  a  partial  pressure  wave  form  is  another  indication  of  the 
versatility  of  the  CC  rotor  concept. 

The  success  of  feedback,  particularly  teeter  rate  feedback,  has  implications 
beyond  the  achievement  of  dynamic  equilibrium  for  simulation  purposes.  The  dynamic 
equilibrium  which  was  reached  with  feedback  was  not  the  same  as  the  dynamic 
equilibrium  reached  through  trial  and  error  without  feedback;  hence,  its  utility  is 
limited  for  the  purpose  of  accelerating  convergence  to  an  open  loop  condition. 
However,  teeter  rate  feedback  was  so  successful  that  it  should  be  considered  as  the. 
primary  means  of  controlling  this  type  of  rotor.  For  example,  rather  than  connect 
a  pilot's  control  stick  mechanically  to  a  pneumatic  valving  system  to  generate  the 
cyclic  pressure  variation  needed  to  control  the  rotor's  thrust  vector  orientation 
(mean  tip  path  plane),  the  pilot's  control  stick  orientation  could  be  used  as  the 
commanded  rotor  tilt  angle  with  the  feedback  control  system  generating  the  appro¬ 
priate  pressure  variation  to  achieve  the  desired  mean  rotor  tilt  angles.  This 
method  would  eliminate  the  need  for  a  mechanical  linkage  between  the  pilot's  control 
stick  and  the  rotor  or  control  system  without  introducing  additional  complexity 
into  the  helicopters  stability  augmentation  system  (SAS) .  (The  SAS  must  modify 
mechanical  control  inputs  to  account  for  changes  in  rotor  phase  relationships  which 
functions  of  speed,  thrust,  etc.)  The  one  complexity  which  would  be  introduced 
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is  the  need  for  stick  position,  rotor  position,  and  rate  sensors.  This  method  of 
controlling  a  rigid  gimballed  rotor  is  very  promising  and  may  have  some  application 
for  articulated  CC  rotors. 

STABILITY 

As  previously  discussed,  the  tests,  applied  to  determine  stability  character¬ 
istics  once  dynamic  equilibrium  was  reached,  consisted  of  monitoring  the  rotor's 
response  to  pu.lse  and  step  gusts.  The  gusts  applied  were  in  the  direction  of  the 
forward  flight  velocity  and  added  either  10%  or  50%  to  the  magnitude  of  the  flight 
velocity.  The  initial  conditions  for  all  cases  were  those  of  open  loop  dynamic 
equilibrium  shown  in  Figure  8.  These  were  50  ft/sec  (15.24  m/s)  forward  flight 
speed,  -4.6  deg  (-0.0803  rad)  blade  collective  angle,  1.8  blade  duct  pressure  ratio, 
and  20%  chord  principal  axis  location.  For  the  open  loop  cases,  the  sine  cyclic 
pressure  ratio  was  -0.0327  and  cosine  cyclic  pressure  ratio  was  0.0384.  One  set  of 
flap  rate  feedback  cases  was  run  with  a  gain  of  k^  =  400  and  with  zero  sine  and 
cosine  cyclic  pressure  ratios. 

Time  history  rotor  dynamics  were  simulated  for  14  rotor  revolutions,  or  about 
3.15  seconds.  The  gusts  were  initiated  at  the  beginning  of  the  second  half 
revolution.  The  pulsed  gust  had  a  duration  of  one-half  revolution. 

Open  Loop  Stability 

The  results  of  the  open  loop  stability  simulation  are  shown  in  Figure  16  for 
no  disturbance,  10%  and  50%  pulsed  gusts,  and  10%  and  50%  step  gusts.  The  responses 
to  the  pulsed  gusts  were  relatively  stable  in  the  sense  that  the  excursions  damped 
out  fairly  quickly;  however,  the  rotor  did  not  immediately  return  to  its  initial 
conditions.  When  the  rotor  was  disturbed  with  step  gusts,  it  departed  aft  and  to 
the  left  of  the  equilibrium  condition.  The  departure  angle  was  approximately 
constant  at  33  deg  (0.58  rad),  measured  from  the  gust  axis,  and  the  departure  rate 
was  proportional  to  the  magnitude  of  the  gust.  Another  feature  of  the  open  loop 
step  gust  response  shown  in  the  figure  is  the  increasing  spacing  of  the  circular 
oscillations  which  indicates  that  the  rate  of  departure  is  increaoing.  From  these 
results,  it  is  evident  that  the  two-bladed  rigid  gimballed  rotors  are  unstable  in 
forward  flight  without  feedback. 
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Closed  Loop  Stability 

The  results  of  the  closed  loop  simulations  with  teeter  rate  feedback,  k^  ■  400, 
are  shown  in  Figure  17  for  no  disturbances,  a  50%  pulse  gust,  and  a  50%  step  gust. 

A  significant  difference  between  the  open  loop  and  closed  loop  simulations  was  the 
presence  of  transitional  dynamics  which  may  be  seen  by  comparing  Figures  16a  and  17a. 
These  dynamics  are  present  in  the  closed  loop  cases  because  of  the  selection  of 
open  loop  dynamics  as  initial  conditions.  The  transitional  dynamics  damped  very 
quickly  as  the  closed  loop  equilibrium  conditions  were  approached. 

The  response  to  the  50%  pulse  gust  was  a  brief  set  of  well  damped  two  per 
revolution  oscillations  followed  by  a  rapid  convergence  to  the  closed  loop  equilib¬ 
rium  dynamics.  The  response  to  the  50%  step  gust  was  well  controlled,  also.  It 
consisted  of  brief  oscillatory  excursions  followed  by  convergence  to  a  new  equilib¬ 
rium  dynamic  condition  pitched  slightly  forward  (nose  down)  from  the  original,  no 
disturbance,  feedback  equilibrium  condition. 

The  step  gust  response  which  was  expected  was  a  shift  to  an  equilibrium 
condition  pitched  slightly  aft  from  the  original  feedback  equilibrium  condition.  A 
slight  speed  instability  is  indicated  by  the  nose  down  pitch  which  would  increase 
the  forward  thrust  component  resulting  in  a  further  increase  in  forward  speed. 
However,  the  fact  that  the  dynamics  damped  to  a  new  equilibrium  quickly  and  that 
the  equilibrium  was  close  to  the  undisturbed  feedback  case,  indicates  a  high 
probability  of  maintaining  control  of  the  vehicle's  motion. 

CONCLUSIONS 

HOVER 

1.  Tip  jet  thrust  may  be  minimized  for  a  given  rotor  thrust  by  selecting  the 
proper  combination  of  rotor  collective  pressure  ratio  and  collective  angle. 

2.  Rotor  dynamic  equilibrium  was  insensitive  to  the  distribution  of  lift  be¬ 
tween  angle  of  attack  and  circulation  control. 

FORWARD  FLIGHT 

1.  Rotor  dynamic  equilibrium  was  attainable  through  the  proper  selection  of 
cyclic  blowing  pressure  ratios. 


21 


2.  The  two-bladed  rigid  gimballed  rotor  was  unstable  without  feedback  control 
as  evidenced  by  its  inability  to  reach  dynamic  equilibrium  after  being  disturbed  by 
step  gusts. 

3.  The  rotor  was  less  unstable  with  CC  lift  than  with  angle-of-attack  lift. 

A.  There  was  an  optimum  feather  principal  axis  of  inertia  location  which 

minimized  the  rate  of  rotor  divergence  and  wobble  amplitude. 

5.  Frequency  of  the  rotor's  wobble  changed  from  two  to  one  per  revolution  when 
the  feather  principal  axis  was  located  too  far  aft. 

6.  Feedback  control,  especially  teeter  rate  feedback,  stabilized  the  two- 
bladed  rigid  gimballed  rotor. 

7.  The  teeter  rate  feedback  stabilized  the  rotor  regardless  of  the  distribution 
of  lift  between  CC  and  angle  of  attack. 

8.  The  teeter  rate  feedback  stabilized  rotor  had  predictable  response  resulting 
in  equilibrium  dynamics  when  disturbed  by  large  gusts. 
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Figure  4  -  Rotor  Forces  and  Moments 
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Figure  6  -  Rotor  Hovar  Performance  Map 
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Figure  8  -  Equilibrium  Dynamics  in  Forward  Fligh 
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Figure  9  (Continued) 
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Figure  10  -  Effect  of  Principal  Axis  Location  on  Rotor  Motions  for  100-Percent 
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igure  10  (Continued) 
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Figure  10  (Continued 
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Figure  10  (Continued) 
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l  Figure  11  -  Effect  of  Principal  Axis  Location  on  Rotor  Oscillation  Frequency 
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Figure  11  (Continued) 
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Figure  11  (Continued) 


Figure  lid  -  100-Percent  Angle-of -Attack  Lift,  Principal  Axis 
Located  at  40-Percent  Chord 
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Figure  12  -  Effects  of  Teeter  Feedback 
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ROTOR  PITCH  ANGLE.  0  (rad) 
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Figure  12  (Continued) 


Figure  12  (Continued) 


Figure  12  (Continued) 


Figure  12  (Continued) 
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Figure  12  (Continued) 


Figure  12  (Continued) 


Figure  12i  -  Teeter  Angular  Acceleration  Feedback;  Gain, 


Figure  13  (Continued) 
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ROTOR  PITCH  ANGLE,  6  (rad) 


:re  13  ^Continued) 


FORWARD  PITCH  - ►  AFT  PITCH 

ROTOR  PITCH  ANGLE,  d  <rad) 


Figure  13  (Continued) 
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Figure  13  (Continued) 
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13  (Continued) 
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igure  13f  -  Feather  Angular  Rate  Feedback;  Gain,  -  -100 


Figure  13  (Continued) 
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Figure  13  (Continued) 
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Figure  13h  -  Feather  Angular  Acceleration  Feedback;  Gain, 


Figure  13  (Continued) 
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FORWARD  PiTCH  - - 

ROTOR  PITCH  ANGLE, 


Figure  14  (Continued) 


Figure  14  (Continued) 


Figure  14  (Continued) 


AVERAGE  ROTOR  PITCH  ANGLE,  8  (rad) 


-0.10 


BLADE  DUCT  PRESSURE  RATIO 

Figure  15  -  Effect  of  Lift  Type  on  the  Average  Equilibrium  Rotor 
Pitch  and  Roll  Angles  with  Teeter  Rate  Feedback 


Figure  16  (Continued) 


Figure  16b  -  10-Percent  Pulsed  Gust  Disturbances 


Figure  16c  -  50-Percent  Pulsed  Gust  Disturbances 
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Figure  16  (Continued) 


ROTOR  PITCH  ANGLE,  9  (rad) 


Figure  1?  (Continued) 


ROTOR  PITCH  ANGLE,  0  (rad) 


Figure  17  (Continued) 
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TABLE  1  -  ROTOR  MODEL  CHARACTERISTICS 


Rotor  Radius 

22.0  ft 

(6.71  m) 

Blade  Chord 

2.933  ft 

(0.894  m) 

Airfoil  Thickness  Ratio 

Root  Section 

0.231 

Tip  Section 

0.153 

Airfoil  Camber  Ratio 

Root  Section 

0.0625 

Tip  Section 

0.0100 

Blade  Twist* 

-8.63  deg 

(-0.151  rad) 

Rotor  Solidity 

0.0849 

Rotor  Weight** 

2 r 817. 3  lb 

(12,531.72  N) 

Rotor  Moments  of  Inertia*** 

T'pll 

4,750  slug-ft2 

(6,640  kg-m2) 

Xp22 

95  slug-ft2 

(129  kg-m2) 

Jp33 

4,815  slug-ft2 

(6,529  kg-m2) 

Rotor  Design  Tip  Speed 

615  ft/sec 

(187  5  m/s) 

*Each  blade. 

**Includes  hub  weight. 

***Principul  moments  of  inertia.  See  Table  2  for  moments  of 
inertia  in  the  rotor  fixed  reference  frame. 


TABLE  2  -  VARIATION  OF  ROTOR  INERTIAS  WITH  PRINCIPAL  AXIS  LOCATION 


Principal  Axis  Location 
(x/c) 

2 

Rotor  Inertias  -  slug-ft 

(kg-m2) 

hi 

I22 

*12 

0.20 

4734.85  (6419.92) 

110.15  (149.35) 

-265.13  (-359.49) 

0.25 

4/39.47  (6426.18) 

105.53  (143.09) 

-221.17  (-299.88) 

0.30 

4743.25  (6431.31) 

101.75  (137.96) 

-177.08  (-240.10) 

0.40 

4748.31  (6438.17) 

96.69  (131.10) 

-  88.63  (-120.17) 

0.50 

4750.00  (6440.46) 

95.00  (128.81) 

0.0  (0.00) 

Notes: 

1.  The  remaining  moments  and  products  of  inertia  were  constants  as 

2  2 

follows:  *  *p33  “  slug— f t  (6529  kg-m  );  I ^  -  0;  and  1^  ■  0. 

2.  Refer  to  Figure  2  for  the  relationship  of  rotor  principal  axis 
location  to  the  rotor  fixed  reference  frame. 
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